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We study a model for the Nematic-Smectic-A (NA) transition in aerogel, and find that even 
arbitrarily weak quenched disorder (i.e., low aerogel density) destroys translational (smectic) or- 
der. Ignoring elastic anharmonicities, but keeping anharmonic couplings to disorder, leads to the 
prediction that there is no "Bragg glass" phase in this system: it is riddled with dislocation loops 
induced by the quenched disorder. Orientational (nematic) order is destroyed as well, as is the 
thermodynamically sharp NA transition, in agreement with recent experimental results. 
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Recent experiments |yj2|] have focussed on the study of 
liquid crystals in the random environment of an aerogel, 
as new paradigm system in which to investigate the gen- 
eral problem of quenched disorder in condensed matter 
systems. 

In this Letter we develop a theory of such systems. We 
find that in a harmonic elastic approximation, even arbi- 
trarily weak disorder destroys both the smectic-A phase 
and the NA transition, and creates unbound dislocation 
loops. The formalism we use to demonstrate this defect 
unbinding is new, powerful, and potentially applicable to 
a wide variety of candidate "Bragg glass" B systems. 

Once dislocations are present, the phase is best char- 
acterized as a nematic in a random tilt field. However, 
subsequent examination of orientational fluctuations in 
this nematic leads to the conclusion that tilt disorder 
destroys the orientational order of the smectic layers as 
well. 

However, a "nematic Bragg glass" and a thermody- 
namically sharp " nematic glass" transition may occur [B , 
in agreement with recent dynamic light scattering exper- 
iments M that show a dramatic slowing down of director 
fluctuation relaxations in liquid crystals in aerogel below 
a temperature T g near the bulk NI transition. 

Near the NA transition the center-of-mass nematogen 
molecular density p(r) develops strong fluctuations dom- 
inated by Fourier components near the smectic ordering 
wavevector qo = n2ir/d parallel to the nematic director 
n. Defining the local (complex scalar) order parame- 
ter ip(r) which distinguishes the smectic-A from the ne- 
matic phase la] via: p(r) = Re[po + e lqo ' r 'i/'(i")], where 
po is the mean density of the smectic, we take for our 
free energy: F = F dG + F dp + F dn , where F dG [yj, n] is de 



F d = / d d 



Gennes free energy (which includes the Frank free en- 
ergy) m1) an d Fdp an d F dn are the new disorder parts, 
which couple the smectic density and the nematic direc- 
tor, respectively, to the quenched disorder of the aerogel 
matrix, and are combined into F d — F dp + F dn : 



-Mr)(p-po) 2 + U(r)p+(g(r)-n) 



(1) 



where both the quenched random T c (St(r)) and the 
quenched random potential U(r) are proportional to the 
local aerogel density pa{y), and the quenched field g(r) 
is random and short-range correlated in direction, with 
strength proportional to the local aerogel density. 

Combining the "random field" energy (first two terms 
of Eq.|l|) with the relation between the smectic order pa- 
rameter tp and the density p, we obtain: 



Fdpty] 



d d r\ 



8t{v)\yb\ 2 +V{r)ip + V*(r)ip* 



(2) 



where V(r) = U(r)e lq ° z . Note that, despite the long- 
ranged correlations of U(r), which arise due to the frac- 
tal structure of the aerogel Q, V(r) has only short- 
ranged correlations. This is because correlations of V 
near q = are related to those of U near q = qoz and 
the aerogel itself has no particular spatial structure at 
the wavevector of the smectic ordering, qo- Thus, the 
correlations of V(r) are short-ranged, and hence we can 
accurately capture the long distance physics by taking 
the real space correlations to be zero-ranged, and write 
V(r)V*(r') = A v S d (r - r') where A v = Cu(q a z) (over- 
bar denotes quenched disorder average). Expanding in 
small deviations from perfect nematic order no = z, writ- 
ing n(r) = z + <5n(r), to linear order in <5n(r), F dn (the 
last term in Eq.[y) becomes: 



F 



dn 



f d d rh(r) 



5n , 



(3) 



where we have defined a quenched random tilt field 
h(r) =g z (r)g(r). 

Since we expect g(r) to have only short-ranged cor- 
relations (with range of order the orientational persis- 
tence length of the silica fibers), the above correlation 
function of the tilt disorder should also be short-ranged. 



Furthermore, it must be isotropic. These considerations 
lead to the following form for the correlation function: 
hi(r)hj(r') = A h S d (r — r') <5y, which is short-ranged. 

While it is tempting to directly analyse the model 
F = FdG + Fdp + Fdn written in terms of the smectic 
order parameter ip, we will not do so here p|. Our mo- 
tivation for this is that such a direct approach is known 
to incorrectly predict the lower critical dimension in, e.g., 
the random field Ising model (as well as, e.g. completely 
missing the existence of the Kosterlitz-Thouless transi- 
tion). Indeed, it proves to also do so here. Instead, 
we proceed by assuming the existence of smectic order 
and writing ^(r) = \ip \e lq ° u< ' r \ with a uniform ampli- 
tude IV'ol = const, and u(r) the local displacement of 
the smectic layers from the perfect periodic order. Using 
this low-T ansatz in Eqs.|l|-y, and integrating over the 
nematic director fluctuations <5n, results in the replace- 
ment Sn — ► Vj_u, everywhere in F[u,Sn] (the Higg's 
mechanism) . This leads to the elastic free energy of the 
disordered smectic-A phase, valid in the long wavelength 
limit, to quadratic order in gradients of u, and provided 
dislocations are confined, 



F[u] = fd d r 



j(d z u) 2 + ^(V 2 ± u) 2 +h(r)-V ± u 



■\tp \ (v(r)e iqoU ^ + V* {r)e- iqou{r) ) 



(4) 



To compute self-averaging quantities, we employ the 
replica "trick" . After replicating and integrating over 
the disorder [Hi, we obtain Z n — f \du n ]e~ F ^ Ua \ with 



J[du a }e- F ^ 
F[u a }= [d d r\jr f^(Vlu a ) 2 + ^(d z u a )'- 



(5) 



a,P=l 



A h 



|Vj_(m q - Uf})\ 2 -~fcos[q a (u a - up)] 



where 7 = \\J)q\ 2 A v /2. We have studied this model us- 
ing the standard momentum shell renormalization group 
(RG) transformation, generalized to allow for anisotropic 
scaling: r± — r' ± e e , z = z'e ue ||; the results are the RG 
flow equations in 3d: 



dl{t) 

di 
dK(£) 

d£ 
dB{£) 
d£ 
dA h {£) 

d£ 



{2 + lu -ry)7-Ai7 2 
(w - 2)K , 
(2-u)B, 

uA h + A 2l 2 , 



(6) 
(7) 
(8) 
(9) 



where for simplicity we set the UV cutoff A = 1, and 
77 = q 2 /(4ir^KB), A 1 = 2n/(3irK), A 2 = <z V(4FQ. 
The statistical symmetry under global rotation requires 
that the disorder generated replica off-diagonal terms to 



be invariant under u a (r) — ► u a (r) + ■ r± . In Eq.|| the 
nonlinearities only depend on the difference between dif- 
ferent replica fields and therefore do not depend on the 
"center of mass" field X^«=i u ai which is therefore a non- 
interacting field. This implies that K and B are not 
renormalized by disorder, |^,0 i.e. their flow equations 
are exact, ignoring (for now) the effects of both anhar- 
monic elastic terms and topological defects loops in u. 
Note that 77, which is simply the Caillc exponent for the 
algebraic decay of smectic correlations in the pure smec- 
tic, is unrenormalized. The recursion relation for the 
proper dimensionless coupling constant 7^ = 27 /(SttK) 
can be easily obtained by combining EqsM|j] 



d£ 



(4 - 77)7 - 



?7 7 2 



(10) 



and, as required, is independent of the arbitrary 
anisotropy rescaling exponent ui. From Eq.nO, we then 
find that for 77 < 4 (large elastic moduli), away from the 
NA transition KB > (7q/2567t 2 , the smectic fixed line 
is unstable to disorder. However, this instability to dis- 
order is stabilized by the nonlinear terms in 7, which 
terminate the flow at a new finite disorder fixed line, 
7* = (4 — T])lr\. This new fixed line then controls a glassy 
smectic-A phase, analogous to the super-rough phase of 
crystal surface on a random substrate MM and the vortex 
glass phase of flux-line vortices in type II superconduc- 
tors §§. 

The flow Eq.fnl also implies that the random-field dis- 
order is irrelevant for 77 > 4. Since the bulk modulus B 
vanishes while K remains finite through Tjvai V diverges 
as T — > TZ A , and hence near the bulk NA transition 
Tata j we are guaranteed to have a range of T over which 
the random-field disorder is irrelevant. However, as we 
will see below, because tilt disorder (A^) is a strongly 
relevant perturbation, the 3d quasi-long-range smectic 
order for 77 > 4 will be converted into short-range cor- 
relations, even when the random field disorder given by 
7 is irrelevant. It is essential to stress that the RG flow 
described above (i.e. relevance for r\ < 4 and irrelevance 
for 77 > 4 of the random-field disorder) survives even de- 
spite the strong relevance and runaway of the random tilt 
coupling A h . 

As can be seen from the recursion relations EqM even 
if the bare A/, = 0, tilt disorder is generated by the 
random-field disorder 7 upon renormalization. In con- 
trast to the 2d random-field XY-model, where the gen- 
erated Ah disorder is only marginally relevant and only 
weakly affects the QLRO order found for A^ = (re- 
placing log phase correlations by log 2 ) H0, for the 3d 
smectic-A phase Ah tilt disorder is strongly relevant. 
The effect of the tilt disorder is controlled by the dimen- 
sionless coupling g = A/,/(i?A 3 ), where A = (K/B) 1 / 2 . 
From the recursion relations EqsJfl-f| we find: dg/d£ = 
2,g + (9?r 3 77 2 /4)7 2 . For 77 > 4, 7(f) — > (as we have 



seen), and so dgjdi = 2g, which is trivially solved to give 
g(£) = g(0)e 2e , where g(0) is a constant. Thus, the tilt 
disorder is strongly relevant. We expect this on physical 
grounds since random tilt disorder explicitely breaks the 
rotational invariance of the smectic-A phase. 

For r\ < 4, 7 — ► 7* > 0. Now, in the 2d random- 
field XY model, the existence of a non-zero 7* in the 
low-T phase implied completely different behavior for 
the tilt coupling g(i) than in the high-T phase. In our 
problem, however, solving EqJq gives g(£) — (g(0) + 
9TT 3 rj 2 ^ 2 /8)e 2i — 9n 3 r) 2 ^ 2 /8, which asymptotically runs 



away to +00 as 



.-,2/' 



in both phases, in exactly the same 



way ! Non-universal constants (like go) change but the 
scaling (e ) does not. This implies that equal time cor- 
relation functions scale in exactly the same way in both 
the "glassy " (77 < 4) and non-glassy (n > 4) phase. We 
will therefore calculate them in the non-glassy (77 > 4) 
phase, where we can set 7 = 0; our results, however, will 
apply to both phases. The only difference between these 
two phases is that the dynamics are slower in the glassy 
phase. 

For 7 = we can calculate anything. The quantity of 

interest C(r±, z) = ((u(r±,z) — u(0, 0)) ) is 



C(r ± ,z)=2 



d 2 q±dq z 
(2tt) 3 



(1 



iq-r\ 



MqMqQ) 

<S rf (q + q') 



(11) 



where the quenched and thermal averaged (u(q)u(q')) = 
(u a (q)u a (q')) = <5 d (q + q')G Qa (q) is expressed in terms 
the replicated correlation function, where no sum on a 
is implied, and the replica propagator G Q/ 3 = 5 a p/T q + 
A h q]/T 2 , with T q = Kq\ + Bq 2 , can be read off from 
Eq.H. This together with Eq.Ogives, in the r — > 00 limit, 

d 2 q±dq z qi (l - e^ r ) 



C = 2A h 

A„ 



32ttBA 3 

r 2 , Ei 



-rl/4A|z| 



(2^) 3 (Kqi + Bq 2 ) 2 

U\\z\e 

-4AL, 
Li 



(12) 



Ei 



\AX\z\ 



2 log 



r± 



where Ei(x) is the exponential integral function, we have 
considered a finite system whose shape is a rectangular 
parallelepiped of linear dimensions L± x Lj_ X L z , L z 
being the length of the system along the ordering (z) di- 
rection, and dropped the subdominant thermal (finite T) 
contribution to C(r±, z). In the usual XL Z « L 2 ^ limit, 
the asymptotic behaviors of C(r±,z) are 



A, 



32ttS 2 A l 



c 



A, 



[4A|z| + r\ log I 



,2_l0g/2VALi 



\\z\ »r'i, 
\L Z « L\ 

A|«|«ri, (13) 
XL Z « L 2 , 



An unusual feature of this result is that even the rel- 
ative displacement of two points with finite separations 



(r±,z) diverge as the system sizes (L±,L Z ) go to infin- 
ity. This is because the mean squared real space orienta- 
tional fluctuations (|5n(r)| 2 ) = (|V^u(r)| 2 ) also diverge 

as L± z — > 00: 



(|Jn(r)| 2 



2A* 



d 2 q±dq z 



<ii 



(2nf (Kqi + Bq 2 Y 



4irB 2 \ 3 



log(min[v / AL z ,Lj 



(14) 



Defining the translational correlation lengths £j_ and 
£ z , as the distances rj_ and z at which C(rj_,z) is of 
order a 2 , where a is a lattice constant, gives: £ z = 
a 2 8irKB/A h , and 



U=4a 



ttB 2 \ 3 



A h log(2^XL z /^) 



1/2 



(15) 



Furthermore, because the liquid crystal in aerogel lacks 
long-ranged orientational order, as well, to obtain the X- 
ray scattering from aerogel, one must powder average; 
the broad ring of X-ray scattering that results has width 
Kpowder = (Cz) _1 = A h /(8nBKa 2 ). 

The orientational correlation lengths £° z 1 can simi- 
larly be defined as the values of L±^ z beyond which the 
mean squared orientational fluctuations (|(5n| 2 ) of Eq. 



47rB- ! \ i /A h 



ae 



^A/2a 2 



get to 0(1). This gives: £[ = ae 

C 2 ° = (a 2 /X)e^ B ' 2x3 / A ^ = (a 2 /X)e^ x / 4a \ Thus, orienta- 
tional order persists out to much larger distances than 
translational order, in the limit of weak disorder where 
all the correlation lengths get large. 

All of the above results apply subject to our two initial 
assumptions, that: (1) dislocations were not generated by 
the disorder, and (2) anharmonic terms in the elastic free 
energy could be neglected. 

We will show now that, if we continue to assume (2) 
(whose validity we will investigate in a future publication 
||), assumption (1) is wrong: in the harmonic elastic ap- 
proximation, dislocations are created even by arbitrarily 
weak disorder. However, they are felt only on length 
scales longer than £° , and hence much longer than the 
translational correlation lengths. Thus, our above calcu- 
lations of these lengths remain valid. 

We can include dislocations in the "tilt only" model, 
i.e. Eq.Hwith V(r) — 0. As discussed earlier, this theory 
correctly reproduces all of the static correlation functions 
in both the glassy and the non-glassy regimes. 

The dislocations are characterized by an integer- valued 
3d vector field m(r) defined on the sites r of a lattice 
connected to the displacement field u via: V x v = m, 
with v = Vw, and a dislocation line continuity constraint 
V ■ m(r) = 0. Standard manipulations Ppl, lead to a 
Coulomb gas theory of these dislocation loop defects: 



a-1 



Kq 



-L D-L 



P i -m i (q)m j (-q) + m(q) • a(-q) 



. (16) 



where the inverse of the pure smectic propagator T q = 

<& + X2 lh P ij (<l) = S ij ~ lt<lf/<lh and a (q) is a Fourier 
transform of the quenched field related to the original 
random phase field h(q) via — ia = qx h/q 2 — (z x q) q- 
hq z (l~\ 2 ql)/{T q q 2 ). 

The partition function for this model is then Z[{h}] = 



Z_/{r 



<r)} ' 



-S[{m}] ) where 5 = Fd / T + £ c /r£ r |m(r)| 2 , 



and the sum is over all integer-valued configurations of 
m, satisfying the dislocation line continuity constraint 
V • m = 0, Hd is given by Eq.ll6], and we have added a 
core energy term E c ^2 r |m(r)| . To proceed, we enforce 
the constraint V • m = by introducing a new auxiliary 
field <fi(r), and introduce a dummy gauge field A to me- 
diate the long-ranged interaction between defects loops 
in the Hamiltonian Eq.[16| obtaining 

Z = II I d<p{r)dA{v) J2 e - s[{m} ^' A] S(V ■ A)S(A Z ) , 



{m(r)} 



S = (3 






m(r) • ( - iV<f>(r) + iA(r) + a(r)) + B c |m| 



+ 



2^ K ql l 



(17) 



Performing the summation of the dislocation loop degrees 
of freedom, replacing the resulting Villain potential by a 
cosine, generalizing to a "soft spin" model described by 
a "disorder" parameter ip — \tp\e 1 ^, we finally obtain 0] 
a complex action S 



*=£ 



(V + iA + a)^*(V - iA - a)t/j + t\tp\ 2 + u|-0| 4 



q 



2Kqi 



A(q)| 2 



(18) 



A complete discussion of the behavior of the above 
model ||] is outside the scope of this Letter; here we are 
only interested in the question of whether the disloca- 
tions loops will or will not unbind at any, even infinites- 
imal amount of disorder. Using replicas and computing 
the disorder-averaged free energy, we find that the lowest 
order contribution to the renormalized dual temperature 
tn comes from the average of the "diamagnetic" term 
SS = Er (d A l 2 > ~W )IVf, which gives 



t R =t + (d-2) 



d d q 

(2^p 



Kq\ A h q 2 z qj 
T, q 2 T 2 



(19) 



where we have used the connection between a(r) and the 
quenched tilt disorder h(r), averaged over h(r) using its 
distribution function, and generalized to d-dimensions. 

The second, disorder, term in this integral dominates 
the first as q — > 0. Indeed, this integral diverges in the 
infra-red for d < 3 



d d q 



2 2 



q 2 (q 2 + X 2 qi) 



4 \2 



Jrf-l 



Q± 



ii 



(20) 



where we used the fact that the dominant regime of the 
integral is q z ~ M±- This divergence implies that t^ is 
driven to — oo (note the minus sign) by the disorder in 
d='i. Indeed, we find in 3d t R = t -Q log(L/a)xO(l), 
where L is an IR cutoff (e.g., the lateral extent of the 
smectic layers) and a is the UV cutoff (e.g., ~ 10A, the 
size of the liquid crystal molecules). 

This implies that the dual (dis-)order parameter is al- 
ways in its ordered phase, which, in turn, implies that the 
dislocation loops of the original smectic model are always 
unbound, thereby destroying the smectic order, even at 
T = 0, for any infinitesimal amount of disorder. 

These results of the harmonic theory imply that for 3d 
disordered smectics the dislocations are unbound, even 
at T = 0, for arbitrarily weak disorder. This means 
that there is no thermodynamically sharp continuous NA 
transition, and the low temperature phase must have a 
finite smectic translational correlation length (even at 
T = 0), once the liquid crystal is put in aerogel. A first 
order transition between nematics with short and long 
(but finite) smectic correlation lengths, or a smeared an- 
alytic crossover, are of course always possible and have 
been observed in experiments of Refs. M and 0], respec- 
tively. We have also shown m that once the disloca- 
tions arc unbound, the static director fluctuations are 
precisely those of a nematic in a random tilt field. This 
implies that the system can be thought of as a nematic in 
such a random field, at all temperatures and all disorder 
strengths. This, in turn, can be shown |3J to imply that 
long-ranged orientational (nematic) order is destroyed as 
well, again in agreement with experiments |lj]. 

Finally, preliminary investigation of anharmonic elas- 
tic effects suggests that they may prevent dislocation 
unbinding and stabilize orientational order. |3j Further- 
more, our results do not imply that there is no thermody- 
namically sharp phase transition in this system analogous 
to the nematic to isotropic phase transition in the pure 
system. Indeed, preliminary investigation suggests that 
a kind of " nematic Bragg glass" phase may exist in these 
systems. This possibility will be further discussed in a 
future publication |3J. 
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